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LAPLACIANS ON METRIC GRAPHS: EIGENVALUES, RESOLVENTS AND 

SEMIGROUPS 

VADIM KOSTRYKIN AND ROBERT SCHRADER 

ABSTRACT. The main objective of the present work is to study the negative spectrum of 
(differential) Laplace operators on metric graphs as well as their resolvents and associated 
heat semigroups. We prove an upper bound on the number of negative eigenvalues and a 
lower bound on the spectrum of Laplace operators. Also we provide a sufficient condition 
for the associated heat semigroup to be positivity preserving. 



1. Introduction 

As suggested by the title, the main objective of the present work is to study the negative 
spectrum of (differential) Laplace operators on metric graphs as well as their resolvents 
and associated heat semigroups. Basic notions related to these operators are summarized 
in Section[2]below. More complete accounts can be found in 1 10 1, 1 121 . 1131 . 

It is well known (see, e.g., [7|, |19|) that there are deep interrelations between pro- 
perties of heat semigroups and spectral properties of their generators. More importantly, 
heat semigroups generated by Laplace-Beltrami operators on Riemannian manifolds carry 
a large amount of information on the geometry of the underlying manifolds. As for metric 
graphs some results in this direction can be found, e.g., in 1221 . [23 1. Nevertheless, a sys- 
tematic analysis of heat semigroups for Laplace operators on metric graphs is still missing. 
In this work we perform a first step in closing this gap. 

Below we will prove an upper bound on the number of negative eigenvalues (Theorem 
13. 7> and a lower bound on the spectrum of Laplace operators (Theorem l3.10> . In particular, 
the upper bound provides a very simple sufficient condition for the Laplace operator to be 
nonnegative. The lower bound improves a recent result by Kuchment 1131 . 

Concerning resolvents and heat kernels we will establish sufficient conditions for them 
to be positivity preserving. To achieve this we will provide a closed expression for Green's 
function in terms of the boundary conditions defining the corresponding Laplace operator. 
As a consequence for a class of boundary conditions associated with positive maximal 
isotropic subspaces (Definition 14.41 below) we prove positivity of Green's function (see 
Theorems 14.61 and 15. 1> . The proof of Theorem 15.11 utilizes walks on a graph, a concept 
introduced in our paper 11 21 as a main technical tool for solving the inverse scattering 
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problem on metric graphs. By standard arguments these results imply that the associated 
heat semigroup is positivity preserving. 

For earlier work on heat semigroups generated by Laplace operators on metric graphs 
and their application to spectral analysis we refer to |2|, |4|, |8 1, |9|, 1 17 1, 1 18 1, \ 22\, |23 1. 
Green's functions have been studied in [2|, |3|, 1 10 1, |20|. 

There is a well-known approach to prove the simplicity of the lowest eigenvalue based 
on the analysis of heat semigroups (see, e.g., 121h . Although we do not pursue this in 
detail, merely as an illustration, for a class of boundary conditions we prove simplicity of 
the lowest eigenvalue of Laplace operators on metric graphs without internal edges (see 
Pi'oposition l6.4l below). 

In Section0we consider Laplace operators on metric graphs with no internal edges, for 
which the heat kernel can be computed explicitly. 

Acknowledgements. The authors are indebted to S. Fulling, H. Hofer, M. Karowski, 
W. Kirsch, M. Loss, and M. Taylor for useful comments. We thank the anonymous re- 
ferees for careful reading of the manuscript and helpful suggestions. The participation 
of the first author (V.K.) at the conference "Quantum Graphs and Their Applications" 
held in the summer 2005 in Snowbird, Utah, USA has been supported by the Deutsche 
Forschungsgemeinschaft and by the National Science Foundation. 

2. Basic Structures 

In this section we revisit the theory of Laplace operators on a metric graph Q. The 
material presented here is borrowed from our preceding papers 1 1 1 and 1 1 2 1 . 

A finite graph is a 4-tuple Q = (V,X, £, d), where V is a finite set of vertices, X is a 
finite set of internal edges, £ is a finite set of external edges. Elements in X U £ are called 
edges. The map d assigns to each internal edge i £ Ian ordered pair of (possibly equal) 
vertices d(i) :— {v%, V2} and to each external edge e S £ a single vertex v. The vertices 
vi =: d~(i) and v-i —: d + (i) are called the initial and terminal vertex of the internal 
edge i, respectively. The vertex v — 9(e) is the initial vertex of the external edge e. If 
d(t) = {v, v}, that is, d~(i) = d + (i) then i is called a tadpole. A graph is called compact 
if £ = 0, otherwise it is noncompact. 

Throughout the whole work we will assume that the graph Q is connected, that is, for 
any v, v' € V there is an ordered sequence {v\ — v, v%, . . . , v n = v'} such that any two 
successive vertices in this sequence are adjacent. In particular, this implies that any vertex 
of the graph Q has nonzero degree, i.e., for any vertex there is at least one edge with which 
it is incident. 

We will endow the graph with the following metric structure. Any internal edge i E X 
will be associated with an interval [0, cti] with a; > such that the initial vertex of i 
corresponds to x — and the terminal one - to i = a,. Any external edge e S £ will 
be associated with a semiline [0, +00). We call the number a, the length of the internal 
edge i. The set of lengths {<2j}j E x, which will also be treated as an element of M' 1 ', will 
be denoted by a. A compact or noncompact graph Q endowed with a metric structure is 
called a metric graph (G,a). 

Given a finite graph Q = (V,X, £, d) with a metric structure a = {a,i}i<=.z consider the 
Hilbert space 



(2.1) H = H(£,I,a)=H s ®Hz, H £ =Q)H e , Hi = ^H l , 



LAPLACIANS ON METRIC GRAPHS 



3 



where Hj = £ 2 (ij) with 

'[0,Oj] if j€X, 
[0,oo) if je£. 



Let ij be the interior of Ij, that is, Ij = (0, dj) if j e X and 7j = (0, 00) if j e £. 

In the sequel the letters x and y will denote arbitrary elements of the product set 
x 

jeeux 

By T>j with j e £ U X denote the set of all ^ € such that V>j (x) and its derivative 
^ (x) are absolutely continuous and tp'J(x) is square integrable. Let T>j denote the set of 
those elements ipj e T>j which satisfy 

^°) = ° for ie£ and ^(0) = = 

^•(0) = tOT and Vj(°)=V'j(%) = 

Let A be the differential operator 

(2.2) (AV) . (*) = jelUf 
with domain 

2?° = D o c ^ 

jesui 

It is straightforward to verify that A is a closed symmetric operator with deficiency indices 
equal to \£\ + 2|J|. 

We introduce an auxiliary finite-dimensional Hilbert space 

(2.3) JC = JC(£, 1) = JC £ © fC { ~ ] © /C^ +) 

with /C f S Cl £ l and /C^ ±) = C' 1 !. Let d JC denote the "double" of JC, that is, d JC=JC®JC. 
For any ^ e I? := (J) Dj we set 

jesui 

(2.4) e d JC, 
with r/> and defined by 

{€(0)}ees 
.{-^i( a i)}i€Xy 

Let J be the canonical symplectic matrix on d /C, 

-(-°, I 

with I being the identity operator on JC. Consider the non-degenerate Hermitian symplectic 
form 

(2.7) W ([0],[V>]):=([0UM>, 

where (-, •) denotes the scalar product in d JC ^ £^\£\+m\), 

A linear subspace M. of d JC is called isotropic if the form vanishes on Ai identically. 
An isotropic subspace is called maximal if it is not a proper subspace of a larger isotropic 
subspace. Every maximal isotropic subspace has complex dimension equal to \£\ + 2\1\. 









(2.5) 


±= {^(0)} ie i , 






\{M a t)}iexJ 
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Let A and B be linear maps of K. onto itself. By [A, B) we denote the linear map from 
d JC = K. © K, to K, defined by the relation 

(A,B) ( X i 0X2) := Axi+Bx-2, 

where xi, X2 G £. Set 

(2.8) :=Ker(A,.B). 

Theorem 2.1. A subspace M C d /C /s maximal isotropic if and only if there exist linear 
maps A, B : K. — > K such that M — M (A, B) and 

(i) the map (A, B) : d JC — > JC has maximal rank equal to \£ | + 2\I\, 

(ii) AS 1, w self-adjoint, AB^ = BA\ 

Definition 2.2. The boundary conditions (A, B) and (A' , B') satisfying i2.9\ are equiv- 
alent if the corresponding maximal isotropic subspaces coincide, that is, M(A, B) = 
M(A',B'). 

The boundary conditions (A, B) and (A' , B') satisfying \2.9\ are equivalent if and only 
if there is an invertible map C : K. — ► JC such that A' = CA and B' = CB (see Proposi- 
tion 3.6 in fHl). 

By Lemma 3.3 in 1 12 1, a subspace M. (A, B) C d JC is maximal isotropic if and only if 
(2.10) M(A, B) 1 - = M(B, -A). 

We mention also the equalities 

M(A,B) ± = [Kcr(A,B)] ± =Ran(AS) t , 

M(A,B) = Ran(-B, A) f . 

In the terminology of symplectic geometry (see, e.g., Section 2.3 in |16|) the equalities 
( 12.1 1> have the following interpretation: The matrix (A, B)^ is a (Lagrangian)/rame for the 
maximal isotropic subspace M.{ A, B) and the matrix (-B, Ay is a frame for At (A, B). 

There is an alternative parametrization of maximal isotropic subspaces of d K. by unitary 
transformations in JC (see 1 1 1 1 and Proposition 3.6 in 1 12 1). A subspace M. (A, B) C d fC 
is maximal isotropic if and only if for an arbitrary k e M \ {0} the operator A + ikB is 
invertible and 

(2.12) 6(k; A, B) := -{A + ikS)" 1 (A - ikB) 

is unitary. Moreover, given any k 6 R\ {0} the correspondence between maximal isotropic 
subspaces M. c d /C and unitary operators S(k; A, B) £ \J(\£\ + 2\X\) on K. is one-to-one. 
Therefore, we will use the notation 6(k; M) for S(k; A, B) with M(A, B) = M. 

There is a one-to-one correspondence between all self-adjoint extensions of A and 
maximal isotropic subspaces of d JC (see 1101 . 1121 ). In explicit terms, any self-adjoint 
extension of A is the differential operator defined by ( 12. 2\ with domain 

(2.13) Dom(A) = {V £ V\ [ip] £ M}, 

where M. is a maximal isotropic subspace of d JC. Conversely, any maximal isotropic 
subspace M. of d JC defines through (12. 1 3I > a self-adjoint operator A(Ai,a). If X = 0, 
we will simply write A(A4). In the sequel we will call the operator A(A4,a) a Laplace 
operator on the metric graph (Q,a). From the discussion above it follows immediately 
that any self-adjoint Laplace operator on Ti equals A(A4, a) for some maximal isotropic 
subspace M. Moreover, A(M,a) = A(M',a) if and only if M = Ml '. 
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From Theorem l2.1l it follows that the domain of the Laplace operator A(A4 , a) consists 
of functions ip e T> satisfying the boundary conditions 

(2.14) A^ + B^/ = 0, 

with (A, B) subject to d2.8i and ( 12.91 . Here ip and ip are defined by (12. 5> . 



With respect to the orthogonal decomposition K = Kg Kj- JC^ any element x 
of /C can be represented as a vector 

({Xe}ee. 
{X^hi 
{x^U 

Consider the orthogonal decomposition 

(2.16) £ = 0£„ 

vev 

with C v the linear subspace of dimension deg(v) spanned by those elements (12. 15l > of K, 
which satisfy 

Xe = if e E £ is not incident with the vertex v, 

(2.17) Xi~ if u is not an initial vertex of i el, 

=0 if v is not a terminal vertex of i el. 

Obviously, the subspaces C Vl and £„ 2 are orthogonal if vi ^ v 2. 

Set d £„ := C v (B C v = C 2dce ^ v \ Obviously, each d C v inherits a symplectic structure 
from d K. in a canonical way, such that the orthogonal decomposition 

vev 

holds. 

Definition 2.3. Given the graph Q = Q(V,T, £ , d), boundary conditions (A, B) satisfy- 
ing (12.91 are called local on Q if the maximal isotropic subspace A4(A,B) of K, has an 
orthogonal symplectic decomposition 

(2.18) M(A,B) = Q) M(v), 

vev 

with M.(v) being maximal isotropic subspaces of d C v . 
Otherwise the boundary conditions are called non-local. 

By Proposition 4.2 in 1121 . given the graph Q — Q(V,T,£,d), the boundary conditions 
(A, B) satisfying ( I2.9i are local on Q if and only if there is an invertible map C : JC — > JC 
and linear transformations A(v) and B(v) in C v such that the simultaneous orthogonal 
decompositions 

(2.19) CA = Q)A(v) and CB = B{v) 

vev vev 

are valid. 

Given a graph Q = (V,l,£,d) to any vertex v 6 V we associate the graph Q v = 
({v},T v , £ v , d v ) with the following properties 

(i) lv = 0, 

(ii) d v (e) — v for all e e £ v , 
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(iii) \£ v \ = degg (v), the degree of the vertex v in the graph Q, 

(iv) there is an injective map ^> v : £ v — > £ Ul such that v £ d o ^ v (e) for all e 6 £ v . 

Boundary conditions (A(v), B(v)) on each of the graphs Q v induce local boundary 
conditions (A, B) on the graph Q with 

A = AO) and B = 



Example 2.4 (Standard boundary conditions). Given a graph Q with \y\ > 2 and minimum 
degree not less than two, define the boundary conditions (A(v), B(v)) on Q v for every 
v *E V by the deg(u) x deg(u) matrices 



A(v) 



(I 








-1 
1 










-1 
1 









-1 

0/ 



B(v) 



/0 







\i i i 



0\ 







1 1/ 



Obviously, A(v)B(v)^ = and (A(v), B(v )) has maximal rank. Thus, for every v £ V the 
boundary conditions (A(v), B(v)) define self-adjoint Laplace operators A(A(v), B(v)) 
on L 2 (Q V ). The corresponding unitary matrices J2. 12b are given by 

(2.20) [6(k; A(v), B(v))] e<e , = ^|^y - 5^ e , 

with 5 e e r Kronecker symbol. 

The boundary conditions B(v)) induce standard local boundary conditions (A, B) 

on the graph Q. 



The following result is Corollary 5 in 1131 . 

Lemma 2.5. Boundary conditions (A, B) satisfying J2.9t are equivalent to the boundary 
conditions (A, B) with 



(2.21) 



A = P, 



Kor B 



B = P, 



Kcr B ! 



where Pxcr B is the orthogonal projection in JC onto Ker B, P^ el B := I — Pkci- b its 
complementary projection, and L : K, — > K, the self-adjoint operator given by 



L — (P|RanSl 

Note that A& = L and Kcr L D Ker B. 



x AP t 



(Kcr B) ■ 



In particular, if KerP = {0} such that L = B~ 1 A, then 6(k; A, B) 
Proposition 3. 19 in 1 12 1). 



L + ik 



(cf. 



Corollary 2.6. Assume that det(A — xB) ^ for some x > 0. Then &(ix; A, B) is 
self-adjoint. If L < 0, then 6(ix; A, B) is a contraction for all x > 0. 

Proof. The assumption dot (A — xB) ^ combined with the fact that the subspace Ker B 
reduces the operator L implies that L — x is invertible. By Lemma l2~5l we have 



(2.22) 



&(ix;A,B) = e(ix;A,B) 



= -Pkc, b - Pk- ct b(L - x)-\L + x)P£ et B 
= -2P Kcr b - (L - x) _1 (L + x), 
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Sje lkX3 for j € 8, 



which is self-adjoint. From \2.22\ it follows that 

||6(ix;yi,B)|| = max{l, \\(L - x) _1 (i + 
If L < 0, then by the spectral theorem || (L - x)~ x (L + x) || < 1 for all x > 0. □ 

3. Eigenvalues of Laplace Operators 

If ip = {V'jljexuf: € Tt is an eigenfunction of the operator —A(A, B,a) corresponding 
to the eigenvalue k 2 el \ {0}, Im k > 0, then it is necessarily of the form 

(3.1) ipj(x;k) = 

The vectors s = {s e } ee £ G Ks, a = {ai} ie i e /Cj , and /3 = {Pi}iei G satisfy 
the homogeneous equation 

(3.2) Z(k;A,B,a) I a 
with 

(3.3) Z(k; A, B, a) = AX(k; a) + ikBY (k; o), 
where 

/I \ /I 

(3.4) X(k;a) = I I and F(k;o) = I -I 

\0 e ik - e _ik -/ \0 — e ik - e _ik -/ 

The diagonal |X| x |X| matrices e ±lk - are given by 

(3.5) [e ±ik ^] jfc = 5 jk e ±ika * for j,k e I. 

The converse statement is also true and we have the following result. 

Lemma 3.1. A i/j given by ( 13. H is an eigenfunction of — A(A4,a) corresponding to the 
eigenvalue k 2 £ R\ {0} if and only if J3.2I possesses a nontrivial solution. The multiplicity 
of this eigenvalue is equal to dim Ker Z(k; A, B, a). 

Proof. If k 2 < 0, the claim is obvious, since ip £ Dom(A(A / f, a)). If k 2 > 0, any 
solution of d3.2i satisfies s = (see the proof of Theorem 3.1 in 1 10 1). Therefore, ?/> G 
Dom(A(M,a)). □ 

As proven in H 1 01 Theorem 3.1] the set of zeros of det Z(k; A, B, a) is discrete. 
Denote 

/0 

(3.6) T(k;a) := e y 

\0 e ik ^ 

with respect to the orthogonal decomposition (12. 3> . 

Theorem 3.2. Assume that I ^ 0. Then A = k 2 ^ racn fnaf det(v4 + ik_B) ^ « an 
eigenvalue of — A(Ai(A, B),a) with multiplicity m if and only if I is an eigenvalue of 

6{k;A,B)T(k;a) 

with geometric multiplicity m. 

If I — 0, then A = — x 2 < is an eigenvalue of— A(A4(A, B)) with multiplicity m if 
and only ifO is an eigenvalue of A — xB with geometric multiplicity m. 
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Remark 3.3. In Q(j) it is shown that det(A + ikB) ^ and det(A - ikB) ^ hold for 
all k > 0. 

Proof. Observe that if det(yl + ikB) ^ 0, then 

AX(k;a) +ikBY(k;a) = (A + ikB)R+(k; a) + (A - ikB)R- (k; a) 
(3.7) = (A + ikB)[l+(^ + ikS)- 1 (A-ikB)T(k;a)]i? + (k;a) 

= (A + ikB)[l- &(k;A,B)T(k;a)]R + (k;a), 



where 
(3.8) 








o > 





I 







e 


-ika , 









f 





: ) 




gika 


0/ 



Hence, by Lemma lXTl | a | G K, is a nontrivial solution to d3.2l > if and only if 

(3.9) [I — 6(k; A, B)T(k\ a)] I a 1=0. 

If X = 0, then equation ( 13 .21 simplifies to (A + ikB)s = with s G K.£. □ 

3.1. Positive Laplacians. Let M. C d /C be a maximal isotropic subspace. Let (A, £?) 
be arbitrary boundary conditions satisfying A4(A, B) = M.. Let n+(AB^) (n_(ABT), 
no (AB^ ) , respectively) be the number of positive (negative, zero, respectively) eigenvalues 
of AB'. By Sylvester's Inertia Law 

n±(CAB*tf) = n±(AB^), n^CAB^C^) = n (AB^) 

for any invertible C : K. — ► K. Since for any such C 

M(CA,CB) =M(A,B), 

we may define 

n±{M) := n±(ABt), n (M) := rc (ABt) 

for arbitrary boundary conditions [A, B) satisfying A4(A, B) = A4. 

We mention several properties of numbers n±(M.) and na{M). First, by ( I2.1Q> . 

(3.10) n±{M^)=n T {M), n (M A -) = n (M). 
Second, for any unitary transformation U in JC one has 

(3.11) n ± { d UM) = n±(M), n ( d UM) = n (M), 
where d U denotes the "double" of U, 



(3.12) d U = 



U 
U 



Finally, we have the following obvious inequality 

(3.13) n (M{A,B)) > max{dimKer A, dimKer B}. 
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The numbers n±{M) and no(M) admit the following equivalent characterization in 
terms of the matrices 6(k; A4). 

Proposition 3.4. Let M. C d JC be maximal isotropic. For any k > the number of 
eigenvalues o/S(k; M.)) lying in the open upper (lower, respectively) complex half-plane 
is equal to n + (M ) ( n_ (M. ), respectively). The number of real eigenvalues of & (k; M. ) is 
equal to no(-M). 

Proof. Consider 

withe = &(k;M). Note that 6(k;M) = &(k;A 6 ,B 6 ). Observing that 

(3.14) A 6 Bl = i-Im 6 = ij-(6 - 
we obtain that 

n±(Im&(k;M)) = n±(M), n (Im6(k;M)) = n (M) 
for any k > 0. □ 

Let M. C d K be a maximal isotropic subspace. For arbitrary ip,i/j S DomA(M,a) 
integrating by parts we obtain 

(3.15) (cp,~A(M,a)i;) n = ]T (<p',i>') Hj + ([<p], Q[ip]W, 

jesux 

where Q = ^ with respect to the orthogonal decomposition d K, = K, K.. 

Using ( 12. 1U it is an elementary exercise to check that the orthogonal projection in d JC 
onto the subspace M. is given by 

Pm = (~jf ) (AAi+BBi)-\-B,A) 

' B\Ajtf + BB^y 1 B -^(AAt +BSt)" 1 A 
-At(i4At +BB'*)- 1 B Ai(AA* +BBi)~ 1 A 

where the block matrix notation is used with respect to the orthogonal decomposition d K. = 
K, © K,. Since [<p] and [tp] in ( I3.15l > belong to M., we obtain 

{W\,QW\) d K = ([<p],QmW\)*k,> 

where Qm '■— PmQPm- Using J3 . 1 6I > we get that 

(3.17) Q M = - (A J 4 t + BB^AB^AA^ + BB^)~ l {-B, A). 

Thus, we have proven the following result. 

Proposition 3.5. For any maximal isotropic subspace M. C d fC 

(3.18) (^-A(M 1 a)^) H = (^^'H + d^QMbP})^, 

jesux 

holds for all ip, if) G Dom A(M , a). 
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Remark 3.6. Observe that Qm depends on the maximal isotropic subspace M. only and 
not on the particular choice of the matrices A and B parametrizing M.. Thus, using the 
parametrization of the maximal isotropic subspace M. by matrices (A, B) defined in (12 . 2 1 i 
we obtain 

(3.19) Qm = - (p^fr) L{1 + L 2 )- 2 (-P£ elB , P Kei B + L) , 

\-T Ker B + L* J 

where we have used the equality 

(3.20) AB f = L. 

Since Ker L D Ker B, the equality Aip + Btp' = implies that L 2 ip + Lifj' = 0. Therefore, 
for any tp S Dom(A(.M, a)) we obtain 

(-•PiCer B > P Kcr B + L) [tp] = -P£, B ± + (P K cr B + ^Mj_ 

= -(I + L 2 )P£ lB l + P Kcl Bf- 

Hence, 

L(I + L 2 )- 2 (-P£ clB , P KclB +L) [0] = -{I + L 2 }- 1 ^. 
Therefore, using ( 13 . 1 91 we obtain that 

holds for all tp, tp G Dom(A(.M, a)), where 

Q^=({) o) 

with the block-matrix notation with respect to the orthogonal decomposition d JC = JC ® 
JC. Together with (I3. 1 8I > this leads to the representation for the sesquilinear form of the 
operator — A(A / t, a) obtained previously by Kuchment in 1 1 3 Theorem9]. 

From ( I3.17> it follows that n±(Q_w) = n T (A4). Thus, by Sylvester's Inertia Law, 
Pi'oposition l3.5l immediatelv implies that — A(A / J, a) > in the sense of quadratic forms, 
if the maximal isotropic subspace M C d K. satisfies the condition n+(M.) = 0. 

Using variational arguments, we prove a more general result: 

Theorem 3.7. The number of negative eigenvalues of— A(M , a) counting their multiplic- 
ities is not bigger than n+(M.). If X — 0, then — A(A4) has precisely n + (A4) negative 
eigenvalues. 

For the proof we need the following simple lemma. 

Lemma 3.8. Assume that (A,B) satisfies (I2.9i . The function f(x) := dot (A — xB) 
has precisely n + (M(A, B)) positive zeroes and n-{M(A 1 B)) negative zeroes (counting 
multiplicity). If k = is a zero of f, then its multiplicity equals Rank B—n + {M.{A 1 B)) — 
n_(M(A,B))<n (M). 

Proof. By Lemma l2~5l the boundary conditions (A, B) are equivalent to the boundary con- 
ditions (A, B) defined in ( 12. 2\\ . Now assume that a n > is a zero of / with multiplicity 
m > 1. By ( I3.20L we have 

(3.21) det(A-xS) = cdet(A- xB) =cdet(i-x) 

with a nonzero constant c S C. Thus, x is an eigenvalue of L with multiplicity m. By 
(13.201 it is also an eigenvalue of AB'. 

The case of negative eigenvalues can be considered in exactly the same way. 
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Observe that by J3.2H the total number of zeroes of f(x) equals the dimension of 
Ran £?t, that is, the rank of B. Thus, x = is a zero of /(x) with multiplicity 

RsDkB-n + (M(A,B))-n-(M(A,B)) 

< \£\ + 2\I\ - n+{M(A, B)) - n_(M{A, B)) = n (M(A, B)). 

□ 

Proof of Theorem \3. 71 Set for brevity m = n + (A4). Assume that — A(A4, a) has at least 
m + 1 negative eigenvalues. By the Min-Max Principle the (m + l)-st eigenvalue is given 
by 

A m+ i = sup inf (i/j,-A(M ) a)tp) n . 

A/CH V'6Dom(-A(X,a)) 
dimAf— m ||'0||^=1 

In particular, this implies that 

(3.22) A m+1 > inf (V,-A(7W,a)V)w 

t/)£Dom(-A(A(,o)) 

Wll«=l 
«/>_LAf 

for any m-dimensional subspace M of 7i. 

Let A/" m C d K. be the subspace spanned by all eigenvectors of Qm with negative eigen- 
values. Since Ker Qm D A4 the inclusion JV m C M. holds. We claim that there is an 
m-dimensional subspace M m C H such that [ip] G N m f° r an y G ■Mn- Indeed, choose 
an arbitrary basis %W , • • • , m N m . For every fee {1, . . . , m} solve the equation 

{-A(M,a) - l)^ k) = with [7/.( fe )] = x ( fe ). 

Let 6> : [0, oo] — > [0, 1] be an arbitrary infinitely differentiable function with 9(x) = 1 for 

all x G [0, 1] and = for all x € [2,oo). Now set i^ k) (x) = (x) for all i G I 
and 

^W(z) = v4 fe) (z)0(z) 

for all e G £. Obviously, = [?/> (fc) ] and G Dom(-A(M,o)). Since are 

linearly independent, we can choose M m as a subspace in H spanned by ip^- 1 ' , . . . , 

Now setting TV = Af m in ( 13.221 1 and using ( 13.181 1 we obtain that A m +i > 0, which is a 
contradiction. 

The second statement of the theorem follows immediately from Lemma l3~8l □ 
3.2. Eigenvalue zero. Obviously, if A(.M, a)ip = 0, then t/j has to be piecewise linear, 

(3.23) i>i(xi) = oti + PiXi, iel and ip e {x e ) = 0, e G £. 

Proposition 3.9. Ifn + (M) = 0, ?/zen dim Ker A(/Vf, a) < n (/Vf). 

Proof. Since n + (.M) = 0, from (13. 171 it follows that Q M > 0. Therefore, if ip G 
Ker A(7W , a), then by Proposition l3.5l we have ifj' = 0. Hence, -0 € Ker A. Therefore, by 

dim Ker A(M, a) < dim Ker A. 
By ( 13 . 1 3I > the dimension of Ker A does not exceed no(M), which proves the claim. □ 

Solutions of the Laplace equation on infinite periodic metric graphs have been stud- 
ied recently by Kuchment and Pinchover in 1141 . The kernel of the Laplace operator 
—A(A4,a) on compact graphs with standard boundary conditions has been studied by 
Kurasov and Novaszek in 1 15 1. 
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3.3. Lower bounds on the spectrum. Corollary 10 in 1 13 1 in a slightly generalized form 
provides the following lower bound 

-A(M,a) > -4||£+||(|£| + 2|X|)max{2||i + ||,maxar 1 }, 

where L+ is the positive part of the self-adjoint operator L = (-B|R anS t) _1 ^-P(Kerfl) x 
defined in Lemma l2~5l Observe that the r.h.s. of this estimate depends linearly on the total 
degree of the graph Q, that is, on the number 

^2 de &( v ) = \ £ \ + 2|I|- The following 

theorem provides a lower bound on the spectrum of Laplace operator, which is uniform 
with respect to this quantity. 

Theorem 3.10. Assume that I ^ 0. Then 

-A(M,a)>- S (\\L + \\) 2 , 
where s(t) is the unique nonnegative solution of the equation 



(3.24) stanh(y) 



= t, a = mina, 

iei 



Obviously, s(0) = 0, s(t) is increasing in t and satisfies the estimate s(t) > t for all 
t > 0. 

Remark 3.11. Assume that I = and n + (A4) > 0. Then by Theorem \3.2\ and Lemma 
\2.5\ the number A = — ||L + || 2 is the smallest eigenvalue of—A(M). 

Proof of Theorem \3. 1 6>l If L + = the claim follows from Theorem l3.7l bv (13.201 . Thus, 
we may assume that L + ^ 0. For brevity we set xq := Observe that by (13.211 . the 

number hq is the biggest positive solution of the equation Aet(A — xB) = 0. 
For arbitrary x > x by J2.22I we have 

(3.25) \\6(ix; A, B)\\ = ^±^1 > i. 

x — x o 

By Theorem l3.2l if — x\ < is an eigenvalue of — A(M.{A y B),a), then either 

det(vl - x x B) = 

or 

det$-&(b<i;A,B)T(bfi;a)] = 
holds. In the first case x\ = hq. If the second case holds, then necessarily 

||6(ixi;A,B)T(ixi;a)|| > 1, 

and, therefore, 

||6(ixi ;J 4,B)||||T(ixi;a)|| > 1. 

Therefore, 

xi < x 2 :=inf{x> 0| ||6(ix; A, B)\\ ||T(ix; a)|| < 1}. 
From (I3.25i it follows that 

||e(ix;A,S)||||T(ix;a)|| < ^±^e^ 

x x 

for all x > x . Therefore, x 2 is not bigger then the unique positive solution of the equation 

-e~ a ~ = 1. 



x — Xq 

It is straightforward to verify that this solution is given by s(xq) and s(xq) > xq. □ 
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4. The Resolvent 

In this section we will study the resolvent of the Laplace operator on a metric graph 
(Q,a). In particular, we show that the resolvent is an integral operator. The structure of the 
underlying Hilbert space Ti \2.\\ naturally gives rise to the following definition of integral 
operators. 

Definition 4.1. The operator K on the Hilbert space Ti. is called integral operator if for all 

j, f 6 £ U I there are measurable functions Kj t ji(-, •) : Ij x Ij, — ► C with the following 
properties 

(i) Kj.j'{xj, )tpji{-) £ L 1 (I ji) for almost all Xj G Ij, 

(ii) ip = Kip with 



(4-1) <M*i)= E / K iA 

The (\2\ + \£\) x (|X| + \£\) matrix-valued function (x,y) i— > K(x,y) with 

[ K { x ,y)]j,f = K j,j'( x j^f) 
is called the integral kernel of the operator K. 



Below we will use the following shorthand notation for \A.\Y . 

ip(x) = / K(x,y)<p(y)dy. 



Lemma 4.2. For any maximal isotropic subspace M. C d K. the resolvent 

(-A(M;a) - k 2 ) -1 /or k 2 e C \ spec(-A(M; a)) 

is the integral operator with the (\X\ + \£\) x (|X| + \£\) matrix-valued integral kernel 
r M{ x iU\ Imk > 0, admitting the representation 

r M (x,y; k,a) = r {0] (x,y, k) 

(4.2) : 

--^x 1 k)Z(k;A,B 1 a)- 1 (A-ikB)R + (k;a)- 1 <P(y,k) T , 

where Z(k; A, B 1 a) is defined in ( 13 . 3I >. the matrix &(x, k) is given by 



i>(:r, k) 



4>{x, k) 
<Mz,k) <j)-{x,k) 



with diagonal matrices <fi(x, k) = dia,g{e lkx j }j^g, (f>±(x, k) — diag{e ±lk:Ej and 

r ' (x, y, k)] jd , = iS j>f — — , Xj , Vj e Ij . 

If&et(A + ikB) ^ 0, then 

r M (x,y;k,a) =r^\x,y,k) 
(4.3) ; 

+ -<S>(x, k)i2+ (kjfl)-^ - e(k;M)T(k;a)}- 1 e(k;M)R + (k;a)- 1 ^y, k) T , 
where R+(k; a) is defined in ( I3.8I >. If I = 0, this representation simplifies to 
r M (x, y, k) = r<°> (a:, y, k) + ^(*, k)e(k; A4)</»(y, k). 
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The integral kernel tm (x, y ; k, a) is called Green 's function or Green 's matrix. 
An alternative representation for Green's function of the Laplace operator on a graph 
with X = has been given in 1 3 1. 



Proof. Since k 2 is not an eigenvalue of — A(A4;a), by Lemma l3~T1 the inverse in J4.2i 
exists. Define the operator M(k) as the integral operator with kernel (14.21 . 

Let £ 7i be arbitrary. Set ip — M(k)ip. Obviously, ip{x) is bounded and ip E H. To 
prove that (— A(M; a) — k 2 ) -1 = M (k), it suffices to show that 

(i) ip G Dom(A(M;a)), 

(ii) (-A(M;a) - k 2 )V = </?, 

(iii) the symmetry relation 

(4.4) r M (y, x;k,a) = r M (x , y; -k, a)* 1 

holds. 

Proof of (i). Obviously, ip € T>, where T> = T)j and T>j denotes the set of all 

jesux 

fj E Hj such that fj(xj) and its derivative f'j(xj) are absolutely continuous and fj(xj) 
is square integrable. Set for brevity 

G(k) := -Z(k;A,B,ay\A-ikB)R + (k;a)- 1 . 

Assume that ipj E Hj vanishes in a neighborhood of Xj = and, in addition, in a 
neighborhood of Xj — a,j if j E X. Then 

^-y^^{ yj )dy 3 = f e-^-Vi^ 3 (y 3 )d Vj 

holds for all sufficiently small Xj E Ij and 

holds for all Xj E Ij sufficiently close to a,j if j EX. A simple calculation leads to 



± = ^i?+(k; a)" 1 / $(y, kfp(y)dy + -X(k; a)G(k) / <P(y, k) T V {y)dy 



and 



t = ^R+teaT 1 f $(y, ^) T v(y)dy - ~Y(k;a)G(k) J G kf^y)^, 

where ?/> and ?/>' are defined by ( 12.51 . Therefore, 

i 

+ Bxj/ = — {A - ikB)i?+(k; a)" 1 J $(y, k) T <^2/)*/ 

-Z(W;A,B,a)G(k) / k) T ^(y)dy = 0. 



2k 

Thus, we proved that Alp + Bip 1 = for all ip in a dense subset of 7i. Therefore, Aip + 
Bip' = for all tp G H, which proves the claim (i). 

Proof of (ii). Assume that ipj E Hj is continuous on Ij for every j E X U £. Standard 
arguments based on the Fourier transform show that 

~~k \i + k ) / ^ Xi - VA vAVi)dVo = <Pj{xj), r, E I,, JESUX. 
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where Ij denotes the interior of Ij . This implies that 




for all xj € Hence, (— A(yVf;a) — k 2 )tp = tp for all tp in a dense subset of 7i. Since 
M(k) is bounded, the claim follows. 

Proof of ( Hi). It suffices to prove J4.4I for the representation J4.3i . For the representation 
14.21 the symmetry relation J4.4I > will follow by continuity from J4.5i . 

The relation r^ ' (y, x, k) = r(°) (x, y, ~ k)t is obvious. Observe that i?+(k;a)^ = 
R+(— k;a). Similarly, 

T(k;a) t = T(-k;a) and $(x, k) f = -k) T (transpose). 
The proof of the identity 

e(k;7W)[I-T(k;a)e(k;M)]- 1 = [I — &(k;M)T(k;a)]~ 1 &(k;M). 
is elementary and left to the reader. Now combining these relations we obtain 

tm (x, y, -k, a) 1 = r (0) (y, sc, k) 

+ ^$(y>)i ?+ (k;a)- 1 [I-6(-k;^) t T(k;a)]- 1 6(-k;A^) t i? + (k;a)- 1 $(x ! k) T . 
Using the symmetry property 

6(-k;X) t = 6(k;M), 

we obtain relation i4.4i . 

Proof of g3J. Recall that if det(A + ikB) ^ 0, then 

(4.5) Z(k;A,.B,a) = (A + ikB)[I — S(k; A, B)T(k;a)]R + (k;a). 

Thus, {O} follows from JO . □ 



Observe that by Corollarv l2.6l for any maximal isotropic subspace .M C rf /C the operator 
&(ix; A4) is self-adjoint for all x > 0. Moreover, 6(ix;A4) is real if and only if the 
matrices (A, B) defining M. via J2.8t > can be chosen real. Indeed, by J2.22i . the matrix 
&(ix; Jvi) is real if and only if L is. 

If n + (A4) = 0, then again by Corollary 12. 61 6(ix; Ai) is a self-adjoint contraction for 
all >f > 0. Therefore, I + (5(ix; A4) is a nonnegative operator whenever x > 0. 

Definition 4.3. For any square matrix C we write C >p {respectively, C >~ 0) if all entries 
of the matrix C are nonnegative (respectively, positive). We write C\ )p C 2 (respectively, 
Ci >- C 2 ) ifC 1 ~C 2 ^0 (respectively, d - C 2 >- 0). 

Definition 4.4. The maximal isotropic subspace M. C d K, is called positive, if there is 
a xq > such that I + 6(ix; A4) )p for all x > xq. It is called strictly positive, 
I + 6(ix;A4) >~ Ofor all x > xq. It is called locally strictly positive, if the boundary 
conditions defined by M. are local in the sense ofDefinition \2.3\ andI + &(ix: A4(v)) >~ 
for all x > xq and all v € V. Here M(v) denotes the the maximal isotropic subspace 
from the orthogonal decomposition (12. 1 8i . 

Obviously, strictly positive maximal isotropic subspaces are positive. If M. defines local 
boundary conditions on the graph Q, then there are permutation matrices II e U(|£| + 
2|J|) such that IIS(i>f-; M)U~ l is block-diagonal with blocks {&(ix; M(v))} ve v (see 
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discussion in II 21 ). Therefore, since every element of the matrix II is either 1 or 0, if M 
is locally strictly positive, then Ai is positive but not strictly positive. 



Example 4.5. Consider the standard boundary conditions defined in Examvle \2.4\ From 
(I2.20i it follows that 

[I + 6(k;M(v))]e, e '= T ^ rT >0 
deg(w) 

for all e, e' £ £ v . Thus, the standard boundary conditions are strictly positive on any 
graph Q v and locally strictly positive on the graph Q. 

Theorem 4.6. Assume that the maximal isotropic subspace Ai is strictly positive. Then 
the Green function r>((a:, y; in, a) >~ for all sufficiently large x > 0. Moreover, if the 
graph has no internal lines (X = 0), then tm{x, y, ix, a) )p for all sufficiently large 
x > whenever the maximal isotropic subspace Ai is positive. 

The case of locally strictly positive maximal isotropic subspaces will be treated in the 
following section. 

For the proof of Theorem l4.6l we need the following lemma. 

Lemma 4.7. Assume that I ^ 0. If the maximal isotropic subspace M. — A4(A,B) is 
strictly positive, then 

(4.6) I + [I - &(ix; M)T(ix; a)]^ 1 &(ix; M) y 
holds for all sufficiently large x > 0. 

Proof. Let x > be the largest solution of the equation det(A — xB) = 0. By Proposi- 
tion 3.11 in 1 12 1 there are positive numbers x\ > x and C > such that 

\\e(ix-M)\\ < c 

holds for all x > x\. Observing that ||T(ix;a)|| < e~ a>c with a := mmaj for all n > 0, 

we obtain that ||6(ix; M)T(ix; a)\\ < 1/2 for all x> x 2 :— max{xi, a -1 log(2C)}. 
Therefore, 

I + [I - &{ix: M)T(ix; a^&^x; M) 

(4.7) ~ 

= I + &(ix;M) + ^[e(i>f;7W)T(i^;a)] n e(i>f;7W) 

n=l 

converges absolutely for all x > Furthermore, for all x > xi the estimate 



J2[&{ix; M)T(ix; a)]"S(ix; M) 

n=l 

oo 

<£||S(ix;M)ir 1 ||T(ix,a)|r 



n=l 

2 

- l-||6(ix;M)T(ix;a 



e(i^X)|| 2 ||T(ix;a)|| ^ ^ 2 „_„ 



holds. Since all matrix elements of I + ©(ix; Ai) are rational functions in x, which are 
positive for sufficiently large x, we obtain the claim. □ 
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Proof of Theorem \4.6\ We start with the case X — 0. If e ^ e', then 

[r M (x,y,ix)ke> = ^-e-^[6(ix;M)} e ,e'e-^' > 
for all x e € I e and all j/ e / G I e i. Now consider the case e = e'. Noting the inequality 

we obtain that 

AX AX 

= ^-e~^ (1 + [6(bq A4)] e , e ) e"^ > 
Ax 

holds for all sufficiently large x by the assumption I + S^ix; M) ^= 0. 

Using Lemma l4~7l the case 1^0 can be treated in the same way. □ 

5. Local Boundary Conditions: Positivity of Green's Function 

In this section we will prove an extension of Theorem l4.6l to the case of locally strictly 
positive maximal isotropic subspaces. 

Theorem 5.1. Assume that X ^ and the graph Q has no tadpoles. If the maxi- 
mal isotropic subspace M is locally strictly positive, then the Green function satisfies 
tm( x , U\ y 0/or all sufficiently large x > 0. 

The proof of this theorem is more involved than that of Theorem l4.6l Unlike the case 
of (globally) positive maximal isotropic subspaces, for locally strictly positive maximal 
isotropic subspace the inequality 

1+ [I - &(ix; M)T(ix; a)]^ 1 &(ix; M) 

in general need not hold for all large x > 0. 

Before we turn to the proof of Theorem 15. II we introduce some notion and auxiliary 
results. 

5.1. Walks on Graphs. A nontrivial walk w on the graph Q from j G £ U X to j' G £ U X 
is a sequence 

(5.1) {j,V0,jl,Vl, ■ ■ ■ , jn,V n , j'} 

such that 

(i) ii , • • • , in e 

(ii) the vertices vq G V and v n G V satisfy vq G d(j), vq € d(ji), v n € d(j'), and 
v„ G 9 (i„); 

(iii) for any k G {1, . . . , n — 1} the vertex G V satisfies G <9(jfc) and Vk G 
9(fcfc+i); 

(iv) Vk — Vk+i for some k G {0, . . . , n — 1} if and only if is a tadpole. 
If j, j' G £ this definition is equivalent to that given in 1 12 1. 

The number n is the combinatorial length |w| com b and the number 

n 
fe=l 

is the metric length of the walk w. 

A trivial walk on the graph Q from j G £ U X to j' G £ U X is a triple {j, v, j'} such 
that i? G d(j) and t> G d(j'). Otherwise the walk is called nontrivial. In particular, if 
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d(j) — {vo,vi}, then {j,vo, j} and {j,vi,j} are trivial walks, whereas {j, Vo,j,vi,j} 
and {j, vi , j, Vo, j} are nontrivial walks of combinatorial length 1. Both the combinatorial 
and metric length of a trivial walk are zero. 

We will say that the walk ( 15. II leaves the edge j through the vertex vq and enters the 
edge f through the vertex v n . A trivial walk {j, v,j'} leaves j and enters j' through the 
same vertex v. 

A walk w = {j, vo, ji,v±, . . . , j n ,v n ,j'} traverses an internal edge i S T if jk — i for 
some 1 < k < n. It visits the vertex v if Vk — v for some < k < n. The score n(w) 
of a walk w is the set {rii(w)}i £ x with rti(w) > being the number of times the walk w 
traverses the internal edge i £ I. 

We say that the walk is transmitted at the vertex Vk if either Vk = d(e) or Vk = d(e') 
or Vk G d(ik), Vk S d(ik+i), and ik ^ ik+i- We say that a trivial walk from e' to e is 
transmitted at the vertex v = 9(e) = d(e') if e ^ e'. Otherwise the walk is said to be 
reflected. The walk is said to be reflectionless if it is transmitted at any vertex visited by 
this walk. 

Let x, x' be two arbitrary distinct elements of the canonical orthonormal basis in /C, that 
is all components of \ are zero with the exception of one which is equal to 1. Let v be 

(i) the initial vertex of the internal edge j S I if x S IC j , 

(ii) the terminal vertex of the internal edge j £ X if \ € K,^\ 

(iii) the initial vertex of the external edge j S £ if x € JCg . 

Assume that v' is determined by the same rule from \' ■ 

Assume that the maximal isotropic subspace M. C d JC defines local boundary condi- 
tions in the sense of Definition l2.3l It is straightforward to check that for any m e N the 
equality 



( X ,S(k;M)(T(k,a)6(k;M))™x'> = E e ik l w l^(k;w) 



all walks w of 
combinatorial length m 
from j to j ' 



holds, where the weight W(k; w) associated with the walk w = {j, v,j\,v\,..., v m -\, 
3m,v',j'} is given by 

T^(k;w) := l&iiK-Miv))}^ 

■ (j[ [&(S>r,M{vi))]n,n + )j [©C*. M (v'))] jm ,y . 

If 77i = 1 the product in the brackets has to be replaced by 1. For m = and v = v' we 
have a similar representation 

( X ,S(k;M)x') = [&{pc,M{v))] jd ,, 

which corresponds to the trivial walk {j, v, j'}. 

Let W^f ' } , a, a' e {+, — } denote the set of all walks from j to j' leaving j through 
d a {j) and entering j' through d a (j 1 ). Observe that these sets are disjoint. 
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Lemma 14.21 and in particular equation (14. 3> imply the following representation for 
Green's matrix 

[rM{x,y,iH-,^)]j,r = [r (0) (x, y,'m)] jd , 

3,3 

+ c -*(a J -x J ) T^(ix;w)c-^l w l e -^' 
(5.2) wew (+ r' 

wew ( :; +) 

3,3 



which holds for all sufficiently large x > 0. 



wew (t ; 

3,3 



5.2. Proof of Theorem 15.11 Below we will present a proof of Theorem 15. II for the case 
j, f G X. If one or both edges are external, the proof follows the same lines and actually 
is simpler than in the case considered below. 

Case I: Assume that j = j'. Set v :— d~ (j). From J4.3I and ( I5.2t it follows that 

[r M V, ix, & )] jd = e- 2 ^ (1 + [6(ix; A(v), B(v))] Jd + 0{^ 2 ^ )) 
holds for sufficiently large h > 0. Since M. is locally strictly positive, we have 

[rM{x,y,\H,^j\j t j > 

for all large x > 0. 

Case //: Assume now that j ^ j'. For arbitrary a G {-,+} we will write 

_ /-, if <r = +, 
\+, if a = -. 

Lemma 5.2. Assume that the graph Q has no tadpoles. Let w G Wj ^ be a walk with 

the smallest metric length among all walks in Wjjf • Assume that w is not reflectionless. 

Then there is a reflectionless walk w' G VVj ^ U Wj ^ ^ U wj ^■', <T ^ /row j fo j' smc/i 
f/za/ 

w = {j,a CT (j),w'} ^ w'Gwgr , 

(5.3) w = {w',cT' (/),/} (f w'GWg'^, 

w = ]./.<r U: .w'.<r' // w' g 

Proof. Since w G Wj ^■', cr ' is a walk with the smallest metric length among all walks in 

Wj "j? , it is reflected at at least one of the vertices in d(j) and d(j'). With this observation 
the claim is obvious. □ 
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Lemma 5.3. Assume that the walk w( CT,<7 ' G VV^'/" ^ is reflectionless. Set 



w (^') := {ij5 . (j))W ( CT ,o } e vvjr , 



w 
w 



Then 



(5.4) 



(+ -)l 


= |w(---) 


| +aj, 


(-,+)! 


= |w(--) 


1 


(+,+)! 


= |w(--) 


\+ cij + a 



F( Xj ,Vf,x) := e-* a *W(ix;w ( -'- ) )e- , * |w l ~'~V^' 

+ e _ 3< (a,- % -) W /( ix . w (+-)) e -^|w(+'->| e -^ 

+ e-^'W(ix; w(-'+5) e -"l w< ~' +) l e -"K'-v) 
+ e -K(o i -«,) W r^ w (+,+)j e -*|w(+'+'| e -K(a J ,-» 1 /) 

is positive for all sufficiently large x > and all Xj € ij, j/j< G 7j/. 

Proof. There are four different cases according to each choice of (a, a'). It suffices to 
consider one case, since the other three cases may be discussed in the same way. We 
pick the case (cr, a') = (— , — ), so the walk w( _,_ ) is assumed to be reflectionless. By 
construction we have 

w 
w 
w 

Set v := d + (j) and v' := d + (j r ). Obviously, 

+-<-))[e(ix;M(v))} jtj , 
(-•-))[6(ix;A4(«'))W, 
(-•-))[6(ix;A4(«))]^[6(ix;MK))]i'^. 
The inequalities 

[&(ix;M(v))] jtj > -I, [&(ix;M(v'))] jt , r > -1 
hold for all large x > since M. is locally strictly positive. Hence we may write 

F(xj,y f) K) =e~ m 'W(ijr,w(-'-))e-"We- m i'(l + a(x))(l + b(>t)) 

where 

-1 < o(x) = e~^*~ Xi) [&(pt\M(v))] jtj , 

holds for all large x. Now W(ix; w^ - '') > since w^ - ' - ) is reflectionless. Therefore 
F(xj,yji , x) is positive for all large x > 0. □ 

Now we are in the position to complete the proof of Theorem l5.ll 

Assume that w is a shortest walk in wj ^ for some cr, a' € {+, — }. By Lemma IS~2l 

this walk is either reflectionless or there is a reflectionless walk w' € VV^f ^ U wj^-f U 

VV^f such that w can be obtained from w' by one of relations J5.3I >. 





°) 


= W(ix; w 






= W(ix; w 


VK(ix;w (+ ' 4 


} ) 


= W(ix; w' 
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Applying Lemma 1531 to the walk w in the first case and to the walk w' in the second 
case we obtain that the corresponding contribution to the sum in d5.2t is positive. Hence, 
the leading term on the r.h.s. of i5.2\ is positive for all sufficiently large x > 0. 



6. Positivity Preserving Heat Semigroups 

Since — A(A4 ; a) is bounded from below, the heat semigroup cxp{i A(A4 , a)} defined 
by the spectral theorem is a bounded operator. 

Lemma 6.1. The heat semigroup exp{iA(.M, a)} associated with the Laplace operator 
— A(M;a) is the integral operator with (\X\ + \£\) x (|X| + \£\) matrix-valued integral 
kernel pt(x,y; M, a). The integral kernel is bounded for every t > 0, that is, for any t > 
there is a constant Ct > such that the bound 

\\pt{x,y]M,a}]j,j>\ < C t 
holds for all x, y G x and all j, j' G £ U T. Moreover, it is infinitely differentiable 

j££UT 
o 

for all x 

o 

Recall that Ij denote the interior of Ij. 

The integral kernel p t (x,y; M, a) is called the heat kernel associated with the Laplace 
operator A(A4, a). If X — we will write pt(x, y; A4). 
For any p E [1, oo] we set 

D>{Q) := = {lifoeeux I ^ G LP(I 3 )}. 

Proof of Lemma |6~71 Choose k > so large that — x 2 < inf spec(— A(.M, a)). Then 

exp{tA{M,a)} = (A(M,a) ~ x 2 )- 1 {A{M , a) - x 2 ) exp{tA(A^, a)}. 

By the spectral theorem (A(vW, a) — x 2 ) exp{tA(Al, a)} is bounded as a map fromi 2 (t/) 
to itself. By LemmaPkHthe resolvent (-A(M;a) - maps L 2 (<?) into L°°{Q). By a 
general theory of Carleman operators, it follows that the heat semigroup exp{tA(A4 , a)} 
is the integral operator with essentially bounded integral kernel. 
To establish smoothness, we note that by a similar argument 

(6.1) {-A M .a + x 2 ) n exp{tA M>!l }(-A Mi a + x 2 ) n 

is a bounded map from L 2 (Q) into L°°(Q) for all n € N and all x > 0. Thus, (16. 1> is an 
integral operator with essentially bounded integral kernel. This implies that 

)n 
\pt(x,y;M,a)]j,j< 

o o 

is bounded for almost all Xj G Ij and yj/ G Ij> . □ 

The self-adjointness of A(A4,a) implies the following symmetry relation of the heat 
kernel 

(6.2) p t (y, x;M,a) =pt(x, y;M,a)K 

For special graphs and special boundary conditions (see |2|, [4|, [8|, [9} as well as 
SectionQbelow) the integral kernel of the heat semigroup can be computed explicitly. 
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Definition 6.2. ip = {ipj}jexu£ € H is called nonnegative, in symbols ip > 0, ifipj(x) > 
/or Lebesgue almost all x G Jj /or any j G X U £. The semigroup exp{iA(.M, a)} is 
called positivity preserving /fexp{iA(.A/f, a)}-0 > holds for all nonnegative ip £ Tt. 

Theorem 6.3. (i) If M. is strictly positive, then the heat semigroup e tA ( M '3d associated 
with the Laplace operator A(M,a) is positivity preserving for all t > 0. If T = it 
suffices to assume that M is positive. 

( ii) If M. is locally strictly positive and the graph Q has no tadpoles, then the heat 
semigroup e tA ( M associated with the Laplace operator A (Ai , a) is positivity preserving 
for all t > 0. 

Proof of Theorem \6.3\ Using the operator-valued Euler formula (7), 1191 

s-lim (--A(M,a) + l) " = s-lim (-Y (-A(M,a) + -Y" = e tA ^ M ^. 

n— too y n J n^oo \ t J \ t ) 

the claim (i) follows from Theorem l4.6l and the claim (ii) from Theorem l5.ll □ 

We close this section with a simple applications of results of the last two sections to 
negative spectrum of Laplace operators. 

Proposition 6.4. Assume thatX = 0. If M is strictly positive, then the smallest eigenvalue 
A < of— A(M) is simple. 

Proof. We call %p = {ipj}j^£ positive, in symbols %p > 0, if ipj(xj) > for Lebesgue 
almost all Xj G Ij for all j G £. Observe that by Lemma \4~2\ and Theorem 14.61 the 
resolvent (—A(A4) + k 2 )^ 1 is positivity improving for all x > ^/—\, that is, 

(-A(M) + x 2 y 1 ip > 

holds for all ip > 0. Applying Theorem XIII. 44 in 1 21 1 we obtain the claim. □ 



7. Some Examples 

Throughout this section we assume that the connected graph Q has no internal edges, 
Q= ({u},0,£,d)withdeg(«) = \E\ > 1. 

For standard boundary conditions (see Example l2.4> the heat kernel has been calculated 
in GO: 

deg(v) -2 

gt(xe-ye) j — -—g t (x e + y e ), e = e, 

(7.1) \p t (x,y;M)] e , e , = \ 2 de S^ 

r9t(xe+ye>), e^e', 



where 



deg(u)' 



9ti x ) = -7= cx P{~ x2 / At }- 



By Theorem [fO] the heat semigroup associated with the heat kernel pt(x, y) is positivity 
preserving (see Example l4.5> . Alternatively, positivity of the heat kernel can be deduced 
directly from (17. II using the inequality gt(x e + y e ) < gt{x e — y e ) for all x e ,y e > 0. 

Below we will derive an explicit representation for the heat kernel for a class of bound- 
ary conditions which to the best of our knowledge has not been treated before. We start 
with recalling the following well-known result (see, e.g., JB |, |6 1). 
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Lemma 7.1. For t > the function ht(s, A), defined in terms of the complementary error 
function 

erfc(iy) = —j= I e~ u du 



as 

ht(s,X) = -XV^id g t (s) exp j - ^Vi^j j erfc (^/| _ 

satisfies the relations 

d t h t (s-X)~d 2 s h t (s;X) = 7 
\(2g t (s)-h t (s;X))-dMs;X) = 0. 

From the well-known asymptotic expansion of the error function 1 1 Equation 7.2.14] 
we obtain the asymptotics 

4Xt 

(7.4) h t (s; A) = —g t (s)(l + 0(t)), 

which holds for small t > and fixed s > 0. 

Assume that the maximal isotropic subspace Ai C d JC is defined by the self-adjoint 
A = H and B = I via relation PI Then 

ff-ik 



6(k;,M) 



# + ik 
such that 

I + M) = 2x(x - F)" 1 . 

Observe that the maximal isotropic subspace Ai(H, I) is positive if either _ff )p or all 
off-diagonal matrix elements of H are positive. Indeed, for sufficiently large n we have 
the absolutely convergent expansion 

oo , 

\ + &(ik-M) = 2^ f — 

If H this sum is nonnegative. If all off-diagonal matrix elements of H are positive, 
then iJ > cl for some c > 0. Thus, 

I + x^if >- 

holds for all h > c. This implies that I + S(i>r; A1) >- for all sufficiently large h > 0. 
For any e, e' E £ consider 

f^ e '\s) := [ht(s, H)] e>e ', 

where the matrix-valued function ht(s, H) is defined by (17. 2> via the spectral theorem. 

Theorem 7.2. Assume that the maximal isotropic subspace M. C d JC is defined by the 
self-adjoint A = H and B = I via relation (12.8V The heat kernel p t {x,y; AA) of— A(At) 
is given by 

(7.5) \p t (x,y;M)] e ,e> = gt(x e - y e ) 5 e ,e> + gt(x e + y e ) S e ,e> - ft ( x e + Ve')- 
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For small t > and fixed x, y the following asymptotics holds 

\p t (x, y, M)] e>e > = gt{x e - ye) $e,e> + 9t{x e + He) $e,e' 

( 7 - 6 ) 4iiJ ee - . w ^, NN 

+ 7 , ' s 9t (Xe+Ve') 1 + Ot ■ 

Proof. Obviously, the integral kernel I7.5i satisfies the symmetry relation ( 16. 21 . Set for 
brevity 

(7.7) v[ e ' e \x e ,y e ,) := g t (x e - y e ) 5 e , e > + 9t(Xe + Ve) S e>e > ~ ft'" \ x e + He')- 

By Lemma |7~T1 for any e, e' G £ the function p t ^ e ' e '(ir e , j/ e ') solves the heat conduction 
equation 

d t V ( t' e '\x e ,y e ,)=dlvl e ' e '\x e ,y el ). 
Therefore, we need only to check that M .1\ satisfies the boundary conditions J2.14I with 
A = H and B = I. 
From M3\ we get 

d s ht(s; H) + Hh t (s; H) = 2g t (s)H. 

Therefore, 

(7.8) d y jl e > e '\y e ,) + H e , e ,,ff' e '\y e .)\ x - = 2g t {y e ,)H e , e , . 

e"E£ 

Equation j7.7l i implies that 

^t e ' e \xe,y e <)\x e =0 = 2g t (y e )5 e ,e> - ft*'" \ye') 

and 

d Xe P^ e,e ) {x e ,y e >)\ Xe = = -d yc j( e ' e \y e >). 
From d7.8l > it follows that 

^ H^n'p'f ' e \x e ,y e ,)\ Xe=0 + d Xe T' t f' e \x e ,y e >)\ XF=Q 

e"e£ 

= 2g t {y el )H e<e , £ H e . e „ ff ^ {y e ,) - d y J^'\y e ,) = 
holds for all y e i > 0. This proves the equality d7.5i . The asymptotics J7.6I > follows from 



In particular, the asymptotics J7.6i implies that for small t > the heat kernel is non- 
negative whenever H )p 0. 

References 

[1] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, Dover, New York, 1972. 
[2] S. Albeverio, F. Gesztesy, R. H0egh-Krohn, and H. Holden, Solvable Models in Quantum Mechanics, 
Springer, Berlin, 1988. 

[3] S. Albeverio and K. Pankrashkin, A remark on Krein 's resolvent formula and boundary conditions, J. Phys. 

A: Math. Gen. 38 (2005), 4859 - 4864. 
[4] H. W. K. Angad-Gaur, B. Gaveau, and M. Okada, Explicit heat kernel on generalized cones, SIAM J. Math. 

Anal. 25 (1994), 1562 - 1576. 
[5] J. D. Bondurant and S. A. Fulling, The Dirichlet-to-Robin transform, J. Phys. A : Math. Gen. 38 (2005), 

1505 - 1532. 

[6] H. S. Carslaw and J. Jager, Conduction of Heat in Solids, Claredon, Oxford, 1960. 

[7] E. B. Davies, Heat Kernel and Spectral Theory, Cambridge University Press, Cambridge, 1989. 



LAPLACIANS ON METRIC GRAPHS 



25 



[8] B. Gaveau, M. Okada, and T. Okada, Explicit heat kernels on graphs and spectral analysis, in J. E. For- 
naess (Ed.), Several Complex Variables, (Proceedings of the Mittag-Leffler Institute, Stockholm, 1987-88), 
Princeton Math. Notes Vol. 38, Princeton University Press, 1993. pp. 364 - 388. 

[9] B. Gaveau and M. Okada, Differential forms and heat diffusion on one-dimensional singular varieties. Bull. 
Sci. Math., n. Sir., 115 (1991), 61 - 80. 
[10] V. Kostrykin and R. Schrader, Kirchhoff's rule for quantum wires, J. Phys. A: Math. Gen. 32 (1999), 595 - 
630. 

[11] V. Kostrykin and R. Schrader, Kirchhoff's rule for quantum wires II: The inverse problem with possible 

applications to quantum computers, Fortschr. Phys. 48 (2000), 703 - 716. 
[12] V. Kostrykin and R. Schrader, The inverse scattering problem for metric graphs and the traveling salesman 

problem, preprint (2006). 

[13] P. Kuchment, Quantum graphs: I. Some basic structures, Waves Random Media 14 (2004), S107 - S128. 
[14] P. Kuchment and Y. Pinchover, Liouville theorems and spectral edge behavior on abelian coverings of 

compact manifolds, Trans. Amer. Math. Soc. (to appear). 
[15] P. Kurasov and M. Novaszek, Inverse spectral problem for quantum graphs, J. Phys. A: Math. Gen. 38 

(2005), 4901 -4915. 

[16] D. McDuff and D. Salamon, Introduction to Symplectic Topology, Claredon Press, Oxford, 1995. 

[17] S. Nicaise, Spectre des reseaux topologiques finis, Bull. Sci. Math., II. Ser., Ill (1987), 401 -413. 

[18] S. Nicaise, he laplacien sur les resaux deux-dimensionnels polygonaux topologiques, J. Math. Pures Appl., 
IX. Ser., 67 (1988), 93- 113. 

[19] E. M. Ouhabaz, Analysis of Heat Equations on Domains, London Math. Soc. Monographs, Vol. 31, Prince- 
ton, Princeton University Press, 2005. 

[20] Yu. V. Pokornyl and I. G. Karelina, On the Green function of the Dirichlet problem on a graph, Soviet Math. 
Dokl. 43(1991), 732-734. 

[21] M. Reed and B. Simon, Methods of Modern Mathematical Physics. IV: Analysisi of Operators, Academic 
Press, New York, 1978. 

[22] J. -P. Roth, Spectre du laplacien sur un graphe, C. R. Acad. Sci. Paris, Ser. I Math. 296 (1983), 793 - 795. 
[23] J.-P. Roth, he spectre du laplacien sur un graphe, in G. Mokobodzki and D. Pinchon (Eds.), Theorie du 

potentiel, (Proceedings of the Colloque Jacques Deny, Orsay, June 20-23, 1983), Lecture Notes in Math., 

Vol. 1096, Springer, Berlin, 1984. pp. 521 - 539. 

Vadim Kostrykin, Institut fur Mathematik, Technische Universitat Clausthal, Erz- 
strasse 1, D-38678 Clausthal-Zellerfeld, Germany 

E-mail address: kost rykin@math .tu~clausthal.de, kostrykin@t-online.de 

Robert Schrader, Institut fur Theoretische Physik, Freie Universitat Berlin, Arnim- 
allee 14, D-14195 Berlin, Germany 

E-mail address: schradergphysik . fu-berlin.de 



